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Abstract 

The regular polyhedra have the highest order of 3D symmetries and are exceptionally at¬ 
tractive templates for (self)-assembly using minimal types of building blocks, from nano-cages 
and virus capsids to large scale constructions like glass domes. However, they only represent a 
small number of possible spherical layouts which can serve as templates for symmetric assembly. 
In this paper, we formalize the necessary and sufficient conditions for symmetric assembly using 
exactly one type of building block. All such assemblies correspond to spherical polyhedra which 
are edge-transitive and face-transitive, but not necessarily vertex-transitive. This describes a 
new class of polyhedra outside of the well-studied Platonic, Archimedean, Catalan and and 
Johnson solids. We show that this new family, dubbed almost-regular polyhedra, can be pa¬ 
rameterized using only two variables and provide an efficient algorithm to generate an infinite 
series of such polyhedra. Additionally, considering the almost-regular polyhedra as templates 
for the assembly of 3D spherical shell structures, we developed an efficient polynomial time shell 
assembly approximation algorithm for an otherwise NP-hard geometric optimization problem. 


1 Introduction 

Regular polyhedra are combinatorial marvels. They are simultaneously isogonal (vertex-transitive), 
isotoxal (edge-transitive), and isohedral (face-transitive). In this context, transitivity means that 
the vertices (or edges or faces) are congruent to each other; and for any pair of vertices (or edges or 
faces), there exists a symmetry preserving transformation of the entire polyhedron which isomet- 
rically maps one to the other. Transitivity plays a vital role is assembly, especially self-assembly. 
For instance, if we consider each face as a building block, then face transitivity indicates that a 
single type of block is sufficient to form a shell-like structure. Similarly, edge-transitivity indicates 
that there is exactly one way to put any two blocks together. 

Viruses, natures smallest organisms, utilize this simplicity by encoding the blueprint (RNA/DNA) 
for a single type of protein. These proteins can attach within copies of the same proteins in specific 
ways. The attaching or binding rules is such that when sufficient concentrations of building blocks 



are available, they combine to assemble into a spherical shell, called a capsid. These spherical shells 
have icosahedral symmetry. However, unlike a icosahedron, which only has 20 faces, more than 
half of the viruses have capsids which are formed by a much larger number of proteins. Caspar 
and Klug [6] first addressed the layout of virus capsids using triangular tiles. A similar class of 
assembly is seen in fullerene like particles, with 12 pentagonal and many hexagonal faces, which 
was first characterized by Goldberg m- 

Recently, several researchers have developed efficient constructions and parameterizations of 
Goldberg-like particles [la CHI uni. For instance, Deng et al. [9] studied extensions of Goldberg’s 
construction to other platonic solids, but their study is not exhaustive in characterization and 
enumeration of all the possible cases. In another recent work, Schein and Gayed m developed 
a numerical optimization scheme that takes a Goldberg-like polyhedra, which by construction 
does not have planar faces and is not always convex, and produces strictly convex polyhedra 
while preserving the edge-lengths. However, the resulting polyhedra no longer have any face- 
transitivity. Also, even though the edges have the same length, they are not strictly congruent as 
their neighboring faces are different- hence making such polyhedra unsuitable for using as layouts 
for assembly. 

In this paper, we introduce the almost-regular polyhedra, which preserve the global polyhedral 
symmetry while offer a denser packing with local symmetries, thus ensuring that only one type of 
building block still suffices. Topologically, these polyhedra are face-transitive and edge-transitive; 
and has at most two distinct types of vertices (each vertex is transitive to the other vertices of the 
same type). In some cases, which we have characterized, the polyhedra may become non-convex 
and numerical optimization is required to keep edge-lengths equal and make the faces as congruent 
as possible. Note that even the non-convex polyhedra of this family remain spherical (i.e., any ray 
emanating from the centroid of the polyhedra will intersect it exactly once). 

We believe that our results will greatly impact research in several areas including the field of 
nano-materials which can self-assemble to create nano-structures with desirable properties. For 
example, gold nanorods for cancer imaging and therapy mm, virus capsids and protein-cages 
for targeted drug delivery [SUEHlESKin]. Advanced scientific computation techniques, to explore 
and automatically predict possible nano-structures that can be formed symmetrically by one type 
of building block (eg. an engineered protein) would accelerate the development of new nano¬ 
shell structures. Our theoretical groundwork would greatly support such extensive computational 
techniques. For instance, we show that using our symmetry characterization, a symmetric assembly 
of n particles can be predicted using a algorithm whose running time is only polynomial in n, even 
though assembly prediction is an NP-hard optimization problem. 

2 Background 

The boundary of a convex polyhedron is homeomorphic to a spherical tiling. While the space 
of all convex polyhedra is not enumerable (uncountably infinite), the sub-classes which exhibit 
one or more types of symmetry and congruency conditions are enumerable. Typical congruency 
conditions considered in this context are face transitivity, edge transitivity and vertex transitivity. 
A polyhedron is face-transitive (isohedral) if all faces of the polyhedron are congruent and are 
transitive. In other words, there exists a symmetry transformation of the entire polyhedron which 
would map any specific face. A, onto another specific face, B. Similarly, it is edge-transitive (isotoxal) 
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if all edges of the polyhedron are congruent and transitive in the same sense as above; and vertex- 
transitive (isogonal) if all vertices of the polyhedron are congruent and transitive in the same sense 
as above. 

Depending on which of the above properties are satisfied, we have the following classes of 
polyhedra- 

• Johnson solid: Each face is a regular polygon. But the polyhedron does not satisfy any of the 
transitivity properties. There are exactly 92 such solids- all of them convex. The subclass of 
Johnson solids which have only equilateral triangular faces are called Deltahedra. 

• Catalan solid: Duals of the Archimedean solids. Catalan solids are convex and face-transitive; 
but not edge-transitive of vertex transitive. 

• Semiregular (also called uniform) polyhedra: These are vertex-transitive and each face is 
a regular polygon (of 2 or more different types). Examples of such polyhedra are the 13 
Archimedean solids, and infinite series of convex prisms and anti-prisms with regular polygons 
as the two parallel faces. 

• Quasiregular polyhedra: This class have vertex and edge transitivity. It is easy to show 
that they can have at most two types of regular faces. There are only two such polyhedra- 
cuboctahedron and icosidodecahedron. 

• Regular polyhedra: A regular polyhedron is face, edge and vertex transitive; and has only 
one type of regular face. There are only 5 regular polyhedra in 3 dimensions. These are the 
Platonic solids: Tetrahedron, Cube, Octahedron, Icosahedron and Dodecahedron. 

These special classes are all enumerated and there are only a finite number of them, except for 
the infinite family of prisms and anti-prisms (which are also enumerable as they map directly to 
the set of natural numbers). 

A note on symmetry A Symmetry group consist of a set of symmetry operations (i.e. trans¬ 
formations which map an object to itself) such that the set is closed under the composition (one 
transformation followed by another) operation. The actions of the polyhedral symmetry groups 
can be expressed as pure rotations around different axes through the center of the polyhedron, 
which we assume to be at the origin without loss of generality. We generally refer to the points of 
intersection of these axes with the polyhedron as the locations of symmetry, and refer to the axes 
as the symmetry axes. For instance, the octahedron have 6 axes of 4-fold rotational symmetr}|^ 
going through the four vertices, 8 axes of 3-fold rotational symmetry going through the centers of 
the faces, and 12 axes of 2-fold rotational symmetry going through the centers of the edges. 

2.1 The Almost-Regular Polyhedra and Their Duals 

We define the family of almost-regular polyhedra as all polyhedra which have global polyhedral 
symmetry, have congruent regular faces, and is face and edge transitive. The first condition means 

^n-fold rotational symmetry, also referred to as the symmetry order n, means that a rotation by 27 r/n maps the 
polyhedron to itself 
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that any such polyhedron must have exactly the same number of rotational symmetry axes with 
the same symmetry orders as any specific regular polyhedron. We refer to these as the global 
symmetry axes and locations. We shall refer to these global symmetry axes as gv-symmetry axes, ge- 
symmetry axes and gf-symmetry axes respectively for axes of symmetry going through, respectively, 
the vertices, edge-centers and face-centers of the regular polyhedron. Additionally, due to the 
congruency conditions, all vertices, faces and edges of an almost-regular polyhedron must have 
locations of local symmetry. Local symmetry operations map the vertices, edges, faces immediately 
neighboring the location of local symmetry to themselves, but may or may not map the remaining 
parts of the polyhedron to itself. These axes will be referred to as Iv-, le-, and If-symmetry axes. 



Figure 1: Illustration of the criterion for almost-regular polyhedra. The polyhedra shown 
in (a) and (d) are not almost-regular. In (a), the global symmetries are preserved, but the local 
symmetries are not (for example, it is not locally 2-fold symmetric around the center of DE). In 
(d), local symmetries are intact (note that DEFG and other creased faces are considered a single 
face), but the global 3-fold symmetries are not (for example, around the center of ABD). 

While this class bear some similarity with the Catalan solids, an important distinction is that 
the Catalan solids allow non-symmetric faces (hence, are not isotoxal), as long as all the faces are 
congruent (for example, the solid in Figure [^a) is a Catalan solid, but is not almost-regular. 

Note that, the transitivity properties that make regular polyhedra suitable for assembly, is 
preserved in almost-regular polyhedra, but now the class is richer and more importantly can model 
structures with more than 20 building blocks. 

The dual of an almost-regular polyhedron has regular faces, is isotoxal and isogonal; but may 
or may not be isohedral (they can hae at most two types of faces). The closest known family is 
the semi-regular polyhedra (duals of Catalan solids) which are also isotoxal and isogonal. But the 
semi-regular polyhedra, which includes the 13 Archimedean solids and the family of prisms with 
regular faces, have at least 2 types of faces. 

2.2 Related Prior Work 

The construction scheme, that we propose in the next section, closely follows the one proposed 
by Coldberg m- The family of polyhedra generated by the Coldberg construction rule m are 
fullerene like structures. Fullerene like structures have icosahedral symmetry (symmetry group of 
the icosahedron), and consists of many hexagonal faces and exactly 12 pentagonal faces. The soccer 
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ball is the smallest example of such structure. See Figure for an illustrative description of the 
construction. 


Regular polyhedron 



Unfolded regular polyhderon 
Or 

The 'Net' of the polyhedron 



(b) 



Figure 2: Illustration of the Goldberg construction. The Goldberg construction involves 
unfolding an icosahedron (see (a) and (b)), and then mapping the unfolded icosahedron onto a 2D 
hexagonal lattice scaled and oriented such that all corners of the unfolded icosahedron (its original 
vertices) falls on the centers of some hexagon of the grid (some example scale and orientations (for 
one triangle only) are shown (c)). Finally, the icosahedron is folded back, along with the hexagonal 
grid etched onto its faces. For example, for the scaling and orientation of the red triangle in (c), 
would result in the tiled icosahedron shown in (d). Notice that the new polyhedron has exactly 12 
regular pentagonal faces where the icosahedral vertices originally were, and many regular hexagonal 
faces. 

Caspar and King [6] proposed a similar approach, but using a triangular lattice, instead of a 
hexagonal one, and required that the corners of an edge of the unfolded icosahedron falls on the 
vertices of the lattice (Figure [^. Since, the triangular lattice is simply the dual of the hexagonal 
lattice, the mapping is essentially the same. But the refolded polyhedron now has only regular 
triangular faces. It has 12 vertices where 5 such faces are incident, and many vertices where 6 
faces are incident- the first set are exactly the original vertices of the icosahedron. Notice that this 
polyhedron is exactly the dual of the one constructed using Goldberg’s method (shown in Figure 
[^d), and also overlayed in Figure [^b)). 

Polyhedra produced by Caspar and King’s construction method are almost-regular^ and the ones 
produced by Goldberg’s are duals of almost-regular. But notice that both Goldberg, and Caspar 
and King focused only on the icosahedral case and considered a specific unfolding onto specific 2D 
lattices, and hence only covers a fraction of the possible almost-regular polyhedra and their duals. 
Separately, Pawley m studied other ways of wrapping different polyhedra using different lattices 
from the wallpaper group. However, Pawley did not provide any theoretical characterization of the 
factors related to the possibility or impossibility of such wrappings. We address this issue in the 
following section. 


3 Characterizing All Possible Almost-Regular Polyhedra 

Both Goldberg and Caspar-Klug constructions can be expressed as unfolding a regular polyhedron 
onto a 2D lattice and then refolding it with the lattice etched onto its faces. Pawley’s wrapping 
idea is equivalent. We call these procedures the unfold-etch-refold method. Here, we prove the 
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Figure 3: Illustration of the Caspar-Klug construction. The Caspar-Klug construction 
involves unfolding an icosahedron onto a triangular lattice scaled and oriented such that all corners 
of the unfolded icosahedron (its original vertices) falls on the vertices of the grid (some example 
scale and orientations (for one triangle only) are shown (a)). Then, the icosahedron is folded back, 
along with the grid etched onto its faces. For example, for the scaling and orientation of the yellow 
triangle in (a), would result in the tiled icosahedron shown in (b). 


conditions that must be satisfied to produce almost-regular polyhedra using the unfold-etch-refold 
idea for any regular solid, unfolded in any way, onto any 2D lattice. 

Shepherd’s conjecture 1201 states that all convex polyhedra have a non self-overlapping planar 
unfolding with only edge-cuts. This conjecture is not proved or disproved yet for all possible convex 
polyhedra. However, it is true for the set of special classes we are interested in. Hence, in principle 
it is possible to unfold one such polyhedron and lay it down on a 2D grid, use the grid to draw tiles 
of the unfolded polyhedron, and then fold it back up to get a tiled polyhedron. However, every 
polyhedron actually have many unfolding. For example, icosahedron have 43380 unique unfoldings. 
Caspar and King’s construction produced almost-regular polyhedra using 1 such unfolding, but 
it is not clear whether other unfoldings would also produce similar almost-regular polyhedra, or 
different types of almost-regular polyhedra, or not be almost-regular. To address this question, we 
characterize the relationship of the local and global symmetries of the almost-regular polyhedra, 
and the etched polyhedra (henceforth called tiling) produced using unfold-etch-refold construction. 

First of all, we prove that the lattice onto which the polyhedron is unfolded must be regular. 

Lemma 3.1. The polyhedra generated by an unfold-eteh-refold using any regular polyhedra and 
unfolded in any way, cannot be almost-regular if a non-regular grid/lattice is used. 

Proof. The unfold-etch-refold essentially wraps the lattice/grid over a regular polyhedron. Hence, 
vertices, edges and faces of the regular polyhedron are suppressed, and a new set of vertices, edges 
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and faces appear (all of which belong to the lattice). Now, the local symmetry condition of the 
almost-regular polyhedron requires that every face be symmetric around its center and be congruent 
to each other. This cannot be satisfied if the lattice itself was not regular, or not symmetric around 
some points. □ 

There are exactly 3 regular lattices in 2D- the square lattice, the triangular lattice and the 
hexagonal lattice. The square lattice has 4-fold rotational symmetries at each vertex and face- 
center, the triangular lattice has 6 -fold and 3-fold symmetries at each vertex and face-center; and 
the hexagonal lattice has 3-fold and 6 -fold symmetries at each vertex and face-center. All of them 
have 2 -fold symmetry on the center of each edge. 

Among regular polyhedra, the tetrahedron, the octahedron and the icosahedron have 3-fold 
symmetries at face-centers and respectively 3, 4 and 5-fold symmetries at vertices. The cube has 
4-fold symmetries at vertices and face-centers. The dodecahedron has 3-fold and 5-fold symmetries 
at vertices and face-centers. 

Now we prove another lemma addressing global symmetry conditions at gf-, ge-, gv-symmetry 
axes. 

Lemma 3.2. To satisfy gf-symmetry conditions, all gf-axes must go through a point of the lattice 
that have cn-fold rotational symmetry where n is the order of rotational symmetry around the 
gf-symmetry axes of the regular polyhedron and c is a positive integer. 

Proof. Recall that an almost-regular must have exactly the same number of rotational symmetry 
axes with the same symmetry orders as any specific regular polyhedron, such that there exists 
a rigid body transformation that perfectly aligns these axes to those of the regular polyhedron. 
Since, the unfold-etch-refold is a wrapping, the expected locations and axes of global symmetry 
of the new tiled/etched polyhedra, and the underlying regular polyhedra are already aligned. For 
example, in Figure [^b). A, B and C are locations of 5-fold global symmetry and D is a location 
of 3-fold local symmetry. Let, D be the location of one axis of symmetry going through the face 
of the regular polyhedron and n be its symmetry order. Depending on the chosen unfolding and 
mapping, a particular gf-symmetry axis can have the following three cases- 

• If a gf-symmetry axis goes through a vertex of the tiled polyhedra, then the tiled polyhedra 
can be n-fold symmetric around that axis only if the lattice have cn-fold rotational symmetry 
around its vertices. 

• If a gf-symmetry axis goes through the center of a face of the tiled polyhedra, then the tiled 
polyhedra can be n-fold symmetric around that axis only if the lattice face is cn-regular. 

• If a gf-symmetry does not go through a vertex or a face-center, then the tiled polyhedra can 
not be n-fold rotational symmetric around the axis irrespective of the symmetry of the regular 
grid. 

Hence, the lemma is proved for any unfolding/mapping. □ 

A regular 2 D grid/lattice can be parameterized using two vectors ci, 62 , and a fixed origin O as 
follows. Let, O be a vertex of the lattice and that O has at least 2 neighbors U and V such that 


7 


the vectors O ^ U and O ^ V are not collinear. Then, defining ei = O ^ U and 62 — O ^ 
every point of the lattice can be expressed as linear combinations hei + ke2 where h and k are 
integers. This defines a coordinate system C with O as the origin and 61,62 as the primary axes 
and a co-ordinate (/i, k) representing points on the 2D plane, such that if both h and k are integers 
then, the point lies on the lattice (is a lattice vertex). 



Lattice coordinate systems 

(- 1 , 1 ) ( 0 , 1 ) ( 1 , 1 ) 

. . 

(-1,0) OTei ;i,o) 

(3rl) 


Triangular 


(l'3) 

4 + * 



(- 1 ,- 1 ) 

Square 


Figure 4: Coordinate systems for the 2D regular lattices. In the figures, O is the origin 
and 61 and 62 are two coordinate axes. Any point (x^y) can be reached by vector xei + ye2 from 
the origin. Coordinates of some of the lattice points are shown. 


Lemma 3.3. A lattice described in the C coordinate system is 2-fold rotationally symmetric around 
a point (x^y) for the following cases- 

• Triangular lattice: if and only if both 2x and 2y are integers 

• Square lattice: if and only if both 2x and 2y are integers 

• Hexagonal lattice: if and only if both 2x and 2y are odd integers, or x — 2{y + 1) — 3A: where 
X, y and d are integers. 

Proof. Triangular lattices have 2-fold (or multiples of 2-fold) symmetries around their vertices and 
edge-centers, both of which satisfy the condition that 2x and 2y are integers. Also, no other point 
satisfy these conditions. Square lattices have 2 -fold (or multiples of 2 -fold) symmetries around their 
vertices, face-centers and edge-centers, all of which satisfy the condition that 2x and 2y are integers. 
Also, no other point satisfy these conditions. Hexagonal lattices have 2 -fold (or multiples of 2 -fold) 
symmetries around their face-centers and edge-centers. The edge-centers satisfy the condition that 
both 2x and 2y are odd integers. The face centers are integer solutions of x, y for the family of 
lines X = 2{y -\-l) — 3d where d is an integer. No other point satisfy either of these conditions. □ 

Lemma 3.4. If a face T of a regular polyhedron is mapped to a 2D regular lattice in the following 
ways, then the part of lattice inside T is n-fold rotational symmetric around the center of T, 
where n is the rotational symmetry around the gf-axes of the regular polyhedron, and the set of 
faces intersected by each edge of T is 2-fold rotationally symmetric around the center of the edge 
(location of the ge-axes of the regular polyhedron). 

1. IfT belongs to either a tetrahedron, an octahedron or a icosahedron and is mapped to trian¬ 
gular lattice, such that all corners have integer coordinates. 












2. If T belongs to either a tetrahedron, a oetahedron or a ieosahedron and is mapped to a hexag¬ 
onal lattiee sueh that all earners fall on faee renters of the lattiee (have integer eoordinates 
(x, y) sueh that x = 2{y -\- 1) — 3d where x, y and d are integers). 

3. If T belongs to a eube and is mapped to a square lattiee, sueh that either all the earners fall 
on vertiees of the lattiee, or all the earners fall on faee-eenters of the lattiee. 


Proof. First, we show that such mapping is possible. 

1. Consider the case when T is an equilateral triangular face being placed on a triangular lattice. 
Without loss of generality, we assume that one corner is placed at the origin and another at 
(/i, k) where h and k are integers. Then, it is trivial to show that the third point can be at 
{h + k, —h). 

2. Consider the case when T is an equilateral triangular face being placed on a hexagonal lattice. 
Note that the hexagonal lattice is simply the dual of the triangular lattice, hence a mapping 
that puts corners of T on vertices of the triangular lattice would put corners of T on face- 
centers of the hexagonal lattice. 

3. Consider the case when T is a square being placed on a square lattice. Without loss of 
generality, we assume that one corner is placed at the origin and another at (/i, k) where 
h and k are integers. Then, it is trivial to show that the other points can be placed at 
(h — k,h -\- k) and (—fc, h). 

When an equilateral triangular face T = ABC is placed on a triangular lattice such that the 
corners are at (0,0), (/i, k) and (h + k, —h), the the center O of T is at which is at 

a face center, or at a vertex (if h — k, oi k — —2h). Hence O falls on a location of 3-fold or 
6-fold symmetry. Hence, gf-symmetry is satisfied. The center (x, y) of any edge would satisfy the 
condition that 2x and 2y are integers, and hence falls on a location that have 2-fold (or 6-fold) 
symmetry and satisfies the ge-condition. 

Similar arithmetic can be applied to prove the theorem for the remaining two cases. □ 

We shall refer to the mapping described in Lemma |3.4| a eompatible mapping. 

Lemma 3.5. For any unfolding of a regular polyhedron onto a regular lattiee, if any faee T is 
eompatibly mapped, then all fares are also eompatibly mapped. And, the etehing inside eaeh faee 
are eongruent. 


Proof. For any unfolding, there must be at least another face adjacent to T and shares an edge 
with it. Let that edge be AB. According to Lemma 3.3, for each other point P belonging to T, 
there exists a point P' produced by rotating P around the center of AB by 180 degrees. Also, 
if P was on a vertex (or a face-center), then P' will also be the same. Since only a rigid body 
motion is applied, the new face T' — ABP'... will also be regular and be congruent to T (i.e. it 
is exactly the unfolded face which was neighboring F along AB). Also, for any point X in T, the 
same transformation would map it to a point X' inside F such that the location of X with respect 
to the corners of T is the same as the location of X' with respect to the corners of Hence the 
etching inside T and F are congruent. By propagation of the same argument, all unfolded faces 
have corners at integer coordinates and are congruent, irrespective of the unfolding. □ 
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Lemma 3.6. For any unfolding of a regular polyhedron onto a regular lattiee, if any faee T is 
eompatihly mapped, then the resulting polyhedron will be almost-regular or a dual 


Proof The polyhedron generated by unfold-etch-refold method have local symmetry, due to the 
lattice being regular. 

Lemma |3.4| showed that within the face T, it satisfies global symmetry around the gf-axes. 
Even after the faces are folded back into the complete polyhedron, the congruency of the all the 
face (Lemma |3.5[ ) guarantees that the remaining faces would also map to each other (including the 
etching inside them). 

Lemma [3^ showed that all faces (and edges) are congruent, and also that the etching is 2-fold 
symmetric around the edges of the face T. Hence when folded back, the etchings from a neighboring 
face T' will match up perfectly (intersect the shared edge of T and T' at exactly the same points), 
and in case there are fractions of a lattice-face inside T, exactly its complement will show up on 
the other side (inside T') of the shared edge, thereby all etched-faces are complete. Topologically, 
the etched-faces that cross the face boundaries and the ones that do not, are identical. 


After folding back, the corners of the faces T,T ',... meet at a point. Note that since all faces 
are congruent (Lemma 3.5) and all corners of each face are also symmetric to each other (Lemma 


3.4), gv-symmetry is satisfied. Moreover, when T was mapped such that the corners fell on vertices 


of the lattice, then the gv-axis is surrounded by exactly n congruent and regular etched-faces, where 
n is the rotational symmetry of the gv-axis, and not the rotational symmetry around the lattice- 
vertices. For instance, mapping a tetrahedron onto a triangular lattice will produce 4 vertices 
where 3 edges are incident, and many more (depending on the scale of the mapping) where 6 edges 
are incident. Hence, two types of vertices appear on the polyhedron, the ones coinciding with the 
gv-axis and the ones that do not. On the other hand, if T was mapped such that the corners fell 
on face-centers of the lattice, then the gv-axis is surrounded by exactly a regular polyhedron with 
nc-fold symmetry, where n is the rotational symmetry of the gv-axis and c is an integer. Other 
etched faces will simply depend on the regularity of the lattice. For example, if an icosahedron is 
mapped to a hexagonal grid, then there would be 12 pentagonal faces and many hexagonal faces. 


Hence, the polyhedron is almost-regular if the corners of the faces fell on lattice vertices, or 
dual if the corners fell on face-centers. Tetrahedron mapped to triangular grid’s face-centers is an 
exception, where the dual construction also produce almost-regular polyhedron. □ 


Figure shows a few examples of constructing almost-regular polyhedron and their duals by 
mapping a face of regular polyhedron on regular lattices in a compatible way. Note that the etched- 
faces that cross an edge of T are geometrically not identical to the ones that do not. The ones 
crossing the boundary have a crease inside them, or if they are flattened, they are no longer regular. 
This is addressed in the Section 14.11 

Now we show that any deviation from a compatible mapping results in a violation of some 
global symmetry condition. 

Lemma 3.7. If a regular polyhedron is mapped to a regular lattiee without fulfilling the eompatible 
mapping eonditions, then the resulting polyhedron eannot be almost-regular. 


Proof If some corners of a face T fall on lattice vertices (or face-centers) and other don’t, then the 
etching inside the face cannot be symmetric around the center of the face. Hence the polyhedron 
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Figure 5: Illustration of the constructing almost-regular polyhedron and their duals. 

Top row shows how placing corners of a polyhedral face on vertices of a compatible lattice produces 
an almost-regular polyhedron. The black lines show the original polyhedron, and the red lines show 
the etching/tiling induced by the lattice. The second row shows and example of placing the corners 
at face centers and producing duals of almost-regular polyhedron. Finally, the bottom row shows 
examples of both primal and dual construction using square lattice. 


is not almost-regular. 


If none of the corners of T fall on lattice vertices (or face-centers), and the center of the face is 
not on a lattice vertex or center of a lattice face, then again the polyhedron is not almost-regular 
(by Lemma 3.2). 


Finally, we consider the case where none of the corners of T fall on lattice vertices (or face- 
centers), but the center O of T is on a lattice vertex or a face-center. This would not violate the 
symmetry around the gf-symmetry axes. But we shall show that it would violate the symmetry 
around the ge-symmetry. Since, we are mapping a regular polyhedron onto a compatible lattice, 
there exists a mapping that would place a face T' such that the center falls at O, and the corners 
lie at integer coordinates. We can generate T by scaling and/or rotating T' around O. Now, we 
shall prove that any scaling or rotation of T that moves the vertices off lattice-vertices violate ge- 
symmetry. The proof depends on the geometry of coordinate system and here we shall only show 
if for the cube mapped to the square lattice case. Other cases follow the same pattern of proof. 


Without loss of generality, we assume that the center O is the origin, and we focus on one edge 
AB of the face T with integer coordinates (xi, ^/i) and (x 2 , 2 / 2 )- Hence, the midpoint ( 7 ( ^i +^2 ^ 
falls on a face-center, vertex, or edge-center. After rotating around Ohy 6 degrees, the new location 
of the center will be (x^+X 2 )sinepy,+y 2 )cose y ^ 2-fold 

location if and only if the numerators are integers, which can only happen if 0 is a multiple of 7r/2 
and in that case, T' coincides with T. 
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Now, if T was scaled by s then the new position of the edge-center C would be 
and again the numerators are integers iff s is an integer. And if s is an integer, then the corners of 
T' would also be on lattice vertices or face-centers. 

Hence, scaling or rotating T around its center O will not keep the (7 on a face-center, vertex, or 
edge-center unless the corners of the transformed face T' fall on lattice-vertices. Finally, the fact 
that a regular lattice is not 2-fold symmetric around any point other than the face-centers, vertices, 
or edge-centers completes the proof. □ 


Finally, we conclude our characterization with the following theorem. 

Theorem 3.8. The polyhedron generated by an unfold-eteh-refold is almost-regular if and only 
if a eompatible mapping of a regular polyhedron onto an unfold-eteh-ref old eompatible lattiee is 
performed. 


Proof The proof follows from Lemmas 3.1, 3.2, 3.6 and 3.7 


□ 


3.1 Parametrization 


For the sake of simplicity of presentation, the following discussion is focused solely on mapping 
the icosahedron onto triangular lattices. Other compatible mappings can be discussed in the same 
manner with almost no difference in the theorems/lemmas presented here except for minor changes 
in counting. The choice to focus on the icosahedral case is primarily due to two reasons: first, it 
has the highest level of symmetry among the regular polyhedra which have a compatible mapping, 
and second, it has applications in modeling viruses, fullerenes etc. 

Let £ be a lattice with origin O and axes H and K. Any point in the lattice is expressed using 
coordinates (/i, k) where both h and k are integers. 

Lemma 3.9. Assuming that one eorner A of the faee T of the polyhedron is mapped to the origin 
O of the lattiee (or the nearest faee-eenter for dual eonstruetions). Then speeifying the position of 
another eompatibly plaeed point B{h^k) is sufjieient to parametrize the entire mapping. 


Proof. Since T is regular, there are exactly 2 possible ways to define the other points of T. Lemma 
|3.4| showed that both of these choices will result in a compatible mapping as long as A and B are 
also compatibly mapped. Also by Lemma 3^, these two are congruent. So, any one of them can 
be chosen arbitrarily. □ 


We mention the following lemma whose proof is immediate. 

Lemma 3.10. Topology of any almost-regular polyhedron or its dual ean be parametrized using a 
tuple <V^C^h^k>, where V is a regular polyhedron, C is a lattiee represented using two axes, and 
h and k are integers. 


3.1.1 Combinatorial Results 

We consider the case where V is an icosahedron whose symmetry group will is /, and C is the 
triangular lattice which will be denoted as . Now, we discuss some properties of the lattice. 


12 











Definition 3.11. We define eaeh triangle of the lattiee as a small triangle and use t to denote 
sueh a triangle. Let us define a triple < i^j^ k > where i and j are integers and k G {+, —Let 
the triangle produeed by the interseetions of the lines h — i, k — j and h + k — i+j + 1 (having 
the vertiees (i, j), {i + 1, j) and {i^j + 1)) be denoted Uj^. Similarly, the triangle denoted tij- has 
vertiees (i, j), (i + 1, j — 1) and {i + l,j), and is produeed by the interseetions of the lines h = i + 1, 
k — j and h + k — i + j. 

The proof of the following lemma is immediate from this definition. 

Lemma 3.12. ti-^j-^ki coineide with ti^j^k 2 only if ii = i 2 , ji — j2 o'o^d ki — k 2 . For any 

small triangle in C?, there exists a triple < i,j,k > sueh that tijk represents that small triangle. 

Through etching, the triangular lattice produces a tiling of a face T (which will be called 
a large triangle in this section) of V where each tile is a small triangle. Now we consider some 
properties of this tiling. 

Assuming A is at (0,0), B is {h,k) such that h and k are integers, the tiling produced by 
on T satisfies: 

• The area of T is ^(/i^ + hk + A;^), which is equal to the area of + hk + k‘^ small triangles. 

• In addition to A, B and C, T includes exactly more vertices of C^. Note that any 

vertex that lie on an edge of T is counted as half a vertex. 

• Each edge of T is intersected by at most 2{h + k) — 3 lines of the form h — c, k — d and 
h + k = e, where c, d and e are integers. 

• The number of small triangles intersected by any edge of T is at most 2{h + k — 1). 

The following are some combinatorial properties of the overall tiled polyhedron: 

• There are exactly 20(/i^ + hk + k‘^) small triangles, and the same number of local 3-fold axes. 

• The 12 gf-symmetry axes are surrounded by 5 small triangles. 

• There are exactly W){h? + hk + k‘^ — 1) vertices (not lying on the gf-axes) with 6-fold local 
symmetry. 

Similar properties can easily be derived for other mappings as well. The important point to 
note is that not only the topology, but also the symmetry and combinatorics are also parameterized 
by only h and k. 


4 Constructing all Almost-regular Polyhedra 

In the previous sections we characterized the conditions that must be satisfied by a unfold-eteh- 
refold protocol to produce an almost-regular polyhedron. The characterization immediately lends 
itself to efficient generation of families of such polyhedra whose topology can be parameterized 
using just two variables (discussed below). Furthermore, the symmetry at global and local levels 
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lets us represent the geometry using a minimal set of points. Finally, we show how these properties 
lead to efficient optimization algorithms for constructing 3D shapes with spherical symmetries. 

Note that given a point P with coordinate (i, j) inside T, there exists two other points Q and 
R such that P, Q and R are 3-fold symmetric around the center D of T. The two points Q and R 
have coordinates {h — i—j^k-\-i) and {h + A: + —h — i) respectively. This can be seen by noticing 
that stepping along the H and K axis by i and j units from A(0,0) is symmetric (around D) 
to stepping in —H + K and —H directions by the same units from B{h^k)^ and stepping in —K 
and H — K directions by the same units from C. We can further extend it to triangles and deduce 
the following. 

Lemma 4.1. If A{hi^ki), B{h 2 ^k 2 ) and C{hs^ks) are three points in the HK eoordinate system 
sueh that /ii,/i 25 ^3, ^2, ^3 are integers and ABC is an equilateral triangle whose eentroid is O, 

then the small triangles symmetrie around 

O. 


Now, we define the minimal set of points or non-redundant set of points S such that no two 
points 5^, Sj G S are symmetric to each other around D, and all points in S lie inside or on T. 
Clearly, |S| = , Note that applying operations on S produces all points inside and 

on T. 


Now we are ready to specify a concrete algorithm for computing almost-regular polyhedra (see 


Figure 1^. 


TilingGen('P, C, h, k) 

Constructs an almost-regular polyhedron using compatible mapping of polyhedron V onto lattice such that the 
scaling and combinatorics are specified hy h,k 

1. Assume that the lattice coordinate system is aligned with the Cartesian coordinate system such that the origins 
coincide and one of the axes is aligned to the X axis, and the other lies on the XY plane. 

2. Place one point A at the origin (0, 0) of the lattice, a second point at (h, k). Compute the other corners of the 
face T. Note that we only need to know the number of vertices n of T. 

3. Compute the location of the center D of the face T. 

4. Let Tc be the set of cyclic symmetry operations around D, such that |T(c| = n. 

5. Initialize empty set S 

6. For each lattice point p inside or on T do 

7. Add p to S if and only none of the transformations in Tc applied to p produces a point which is already 
in S. 

8. Compute the transformation Tmap which maps the face T to a face of the polyhedron V. Tmap is composed 
of Tmap^Tmaps'^mapA such that Tmap a translates T along the lattice to take D to the origin O, Tmaps 

a scaling that resizes T to the size of the faces in 7^, then Tmap^ a translation along Z-axis by an amount 
equal to the distance from the center of P to a face-center. 

9. Let Tp be the set of global symmetry operations (from the symmetry group of V. 

10. Define a set of transformations = {T 2 TmapTi\T 2 G Tp&Ti G T§}. 

11. All points on the almost-regular polyhedron is now generated by simply computing Taii{S). 

Figure 6: TilincGen: Algorithm for constructing an almost-regular polyhedron using compatible 
mapping 
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Theorem 4.2. The algorithm TilingGen constructs a minimal geometric representation of the 
almost-regular polyhedron in terms a set of points S embedded onto the XY plane and a set of 3D 
transformations Tair 


Proof Follows from the definition of S and Lemma 3.6 


□ 


Some polyhedron generated by applying TilingGen are shown in Figure 

Note that if the tiles that cross the boundaries of a face T of V are not regular, they would look 
like they have a crease along the edge of the T by definition of the unfold-etch-refold technique. 
Tiles generated by this algorithm will also have the same problem and such tiles will be non-regular, 
and in some cases even non-planar. In the next section we address this issue. 





T= icosahedron 
£ = triangular 
h = 


T= octahedron 


£ = triangular 
h= l,k = 5 


T = icosahedron 
£ = hexagonal 
h = 3, k = 1 


c 


T =icosahed 
£ = triangula 
h = 3,k = 0 


T = tetrahedron 


£ = triangular 
h= l,k=5 



T= icosahedron 


£ = hexagonal 
h = 8,k=8 


Figure 7: Some polyhderon generated by applying TilingGen. 

Some polyhedron generated by applying TilingGen are shown in Figure 

Note that the tiles that cross the boundaries of a face T of V may have a crease along the edge 
of the T. This is a result of the unfold-etch-refold technique. Tiles generated by this algorithm will 
also have the same problem and such tiles will be non-regular, and in some cases even non-planar. 

4.1 Curation of Tiles 

As mentioned before, sometimes the lattice faces, which corresponds to tiles/faces of the generated 
almost-regular polyhedron, crosses the boundary of the polyhedral face T embedded on the lattice. 
During folding, these faces get warped. There can be exactly three types of scenarios (for mapping 
to a square or triangular lattice). 

1. If one corner of T is at (0, 0), and the other corner is at (/i, k) such that either /i = 0 or /c = 0, 
then no lattice face crosses the edges of T, and no curation is required, (see Figure top 
row). 
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2. If one corner of T is at (0,0), and the other corner is at (/i, fc) such that h — k then some 
lattice faces are exactly bisected by the edges of T. The curation is quite trivial in this case. 
li h — the center of T would lie on a lattice vertex. Let the center be and the face T 
be ABC. Then, folding along AD^ BD and CD will not warp any lattice face. Additionally, 
connecting D to the centers of the neighboring polyhedral faces T would satisfy all global 
symmetry conditions as well. This folding will produce a base polytope which actually looks 
like the almost-regular polyhedron with /i = /c = 1. In fact, for any integer i, to polyhedron 
generated for /i = A: = i is nothing but subdivisions of the faces of the /i = A; = 1 polyhedron. 
Interestingly, in some cases, the new folding produces a polyhedron with base geometry like 
some Catalan solids, but will unlike Catalan solids, these will have regular faces and may be 
non-convex. For example, the h = k = 1 polyhedron have the same topology as the pentakis 
dodecahedron (see Figure [^middle row). 

3. If one corner of T is at (0,0), and the other corner is at (/i. A:) such that h ^ k^Szh^k > 0; 
then, for all mappings on the hexagonal lattice, some lattice faces shall cross the edges of 
T in variable ways. There does not exist any folding which can avoid the crossing while 
maintaining global symmetry (as the lines will not meet at the center of the face T). See 
Figure bottom row. In this case, no exact solution exists, and we provide a numerical 
approximation which maximizes the regularity while ensuring that global symmetries are not 
violated (see below). 


One corner is 
assumed to be 
at (0,0) 

Second corner 
is at (h,k) 



or k = 0 
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h = k 


If 

h^k 
h, k>0 


Red lines show the polyhderal face 




































































(0, 

)) 


( 

1,0) 












r " 



N 







/ 



\ 





/ 














/ 

‘h,h 




\ 



/ 





1 


\ 










(0, 

j) 

























































h. 














(0,0 

) 
















No lattice face crosses 
the polyhedral edges 
No curation needed 


Some lattice faces cross polyhedral edge and gets 
warped (symmetrically) 

• Can be curated by folding 
along different edges (dashed 
lines) 

• e.g. Icosahedron is folded 
like an almost-regular 
polyhderon with h=l, k=l 


Lattice faces cross the 
polyhedral edges in 
different ways 
Numerical curation 
required 



Figure 8: Different cases of warping of tiles, and their curations. 

We should mention that recently Mannige and Brooks [1^ explored such pseudo-irregularities 

in icosahedral tilings and suggested the existence of three classes depending on the values of h and 
k. Here, we have generalized that to all almost-regular polyhedra, and quantified the exact number 
of faces that gets warped. 
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4.1.1 Curat ion as a Numerical Optimization Problem 

The goal of curation would be to make all the faces as regular as possible without violating global 
symmetry. A similar problem specifically for the family of polyhedra generated by mapping an 
icosahedron onto a hexagonal lattice was addressed in m- In that work Schein and Gayed aimed 
to make all the hexagonal faces that cross the face boundary, and become creased/non-planar, into 
planar ones while keeping the edge lengths equal. They also showed that it is possible to provide an 
efficient numerical solution to the problem which ensures that no two hexagonal/pentagonal tiles 
lie on the same plane and the overall polyhedron is convex. However, the shapes of the hexagons 
are allowed to get distorted such that the angles are no longer equal, and may vary a lot within the 
same hexagon. Hence, the faces are no longer congruent (or even nearly congruent) to each other. 
This makes such a polyhedron non-amenable for modeling structures formed using a single type of 
building block, for instance viral capsids. In contrast, we want to maintain the congruence of the 
tiles as much as possible. 

When mapping a polyhedron onto a triangular lattice, the generated polyhedra falls under the 
class called deltahedra, polyhedra whose faces are all equilateral triangles. Even though there are an 
infinite family of deltahedra m (our families are also infinite), it has been known since Freudenthal 
and van der Waerden’s work m, that there are exactly eight convex deltahedra, having 4, 6, 8, 
10, 12, 14, 16 and 20 faces; among them only three are regular or have symmetries like the regular 
ones. So, our family of almost-regular polyhedra cannot be convex and regular at the same time. 
We prioritized regularity. 

• Let the set of all points on the generated polyhedron be Sail — ^alli^)- 

• Let El be the set of lattice/tile edges on the generated polyhedron 

• Let E 2 be the set of diagonals of all the tiles on the generated polyhedron 

• Let, for any point p G Saip the functions s{p) and t{p) returns respectively a point g G S and 
a transformation T G Tail such that p — T{q). 

• Let dist{u^ v) be the square of the Euclidean distance between two points. 

In our calculations we shall only update the positions of the points in S, and all other points 
P ^ '^all oil fhe polyhedron will be generated as t{u){s{u)). This ensures that the points are always 
moved in a symmetric way with respect to the global symmetry axes. Hence global symmetry is 
never violated. 

Let be the initial positions of the points in S. As we update the locations of the points in S in 
our algorithm, the squared displacement of each point p G S will be defined as 5{p) = {dist{p^p^)Y^ 
where p^ G is the initial position of p. Also, the squared length of a line segment e{u^v) G 
El n E 2 will be computed as dist{t{u){s{u))^t{v){s{v))Y and be denoted /(e). Let us also define 
Ml = ^EeeEi(Ke)) and /i 2 = ^ EeeE 2 (^(e))- 

Finally, we define an energy function J^(S) as follows: 

^ l^(EeGEi(^(®) “ Ml)) + ]^(EeeE2(^(®) “ M2)) + ]|| (EpeS ^(m))’ 

Now, we minimize the function over the positions of the points in S. 
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Figure 9: Numerical curation of warped faces. Top row: (a) shows a polyhedron with 
icosahedral symmetry and 260 tiles. The color of the triangles are determined by the ratio of the 
longest and the shortest edge of the tile. Ratios 1 to 1.3+ are colored using white to red gradient. 
The triangles at the corners and near the edges of the icosahedron have higher distortion, (b) shows 
the result of numerical curation. All the triangles are now regular, the worst ratio being 1.034, 
and the symmetry is preserved. Bottom row: (c-d) show a similar before-after figure for a smaller 
polyhedra, but one where the warping is more apparent. The numerical optimization brought down 
the ratio of the worst triangle from 2.13 to 1.004. (e) shows a superposition of the two states to 
highlight that points are updated symmetrically. 


18 






This is clearly a quadratic optimization problem over h? + hk + variables, and for most 
practical values of h and k it can be solved efficiently using any numerical solution techniques. We 
chose to use the limited memory variant of Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm 
[4] due to its faster convergence rates. Also since first and second derivatives (hessian) of the energy 
function are straight-forward to compute, the numerical solution does not require finite-differences 
and heuristic based hessians, making the solution under BFGS more efficient and stable. Figure 
shows some examples of numerically curating warped faces. 


5 Constructing Shell Structures 

Viruses, as discussed before, have icosahedral symmetric shells formed by multiple copies of the 
same protein. While several existing works (e.g. [SIEZIIMIIISIEIIIIEIIHI) leverage symmetry to 
analyze a given shell structure, we are the first to propose a generative algorithm to predict all such 
structures. 

We propose that each virus shell is templated on a particular almost-regular polyhedron. We 
consider the assmebly prediction problem where the number of proteins and the structure of a 
individual protein is known, and the structure of the whole shell is unknown and must be predicted. 
We propose the following generation algorithm (Figure [T^ to solve this problem. Note that, while 
the algorithm specifies icosahedral symmetry as that is the relevant one for viruses, it can easily 
be extended to handle other cases. 


ShellGen{n^ B) Given integer n and the structure (and relevant parameters) of building block S, construct a shell 
with icosahedral symmetry. 

1. Find integer solutions for h and k such that c(/i^ hk k‘^) = n, where c is the local symmetry order of the 
face of the tiling (3 for icosahedral). 

2. If no integer solutions are found 

3. Terminate. 

4. Else 

5. Construct a meshed polyhedra or a tiling r with icosahedral symmetry and the solved values of h, k. 

6. Construct a cyclic symmetric super-block SB using c copies of B. 

7. Symmetrically (see Appendix for detail) place a copy of SB on each tile of r. 

8. Search over the space of in-plane rotations and radial translations (symmetrically applied to each SB). 

9. Report the highest scoring model. 

Figure 10: ShellGen: Algorithm for constructing a shell structure template upon an almost- 
regular of requisite size. 

Analysis of the algorithm follows. 

5.1 Datails of Decoration Rules for Constructing Thick Shell Structures 

We assume that the structures of tiles with the same internal symmetry (e.g. small triangle) 
remain the same across the shell. And all tiles are decorated by identical building blocks in a 
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cyclic symmetric configuration. In other words, we assume a direct mapping between the internal 
symmetry of the tile and the set of building blocks used to decorate it. So for example, we decide to 
use a complex of three identical blocks in a configuration (henceforth called a c-tile) to decorate 
a small triangle of the layout, as opposed to using 3 independent blocks to decorate each corner of 
the c-tile. 

Representation of a tile Each tile is represented using five parameters (ui,U 2 ,/). 

• Ui is a unit vector pointing from the origin to the center of symmetry of the tile. 

• U 2 is a unit vector pointing from the center of symmetry of the tile to one representative 
corner of the tile. 

• / is the order of symmetry of the decoration to be placed inside the tile. 

Representation of a c-tile Each c-tile is represented using five parameters (vi,V 2 ,c,o). 

• vi is a unit vector representing the symmetry axis. 

• V 2 is a unit vector orthogonal to the symmetry axis pointing to the centroid of one copy from 
the center of symmetry. 

• c is the center of symmetry. 

• o is the order of symmetry. 

Decorating rules To decorate a tile (ui,U 2 ,/) using a c-tile (vi,V 2 ,c,o), the following must 
be satisfied- 

• o = /. 

• c is translated such that it lies on the line mui — 0 where m is a scalar. 

• Rotated such that vi aligns with ui and V 2 x U 2 is parallel to ui. 

Further, we ensure the following- 

• If Ci and Cj are any two c-tiles of the same order, 

— The c-tiles are translated in a symmetric way. i.e. rrii — rrij. 

— The c-tiles are rotated in a symmetric way. i.e. C 2 ■ • U 2 ^ = V2j * U 2 j . 

The following result is immediate. 

Lemma 5.1. The decoration rules are necessary and sufficient for preserving local symmetry (in¬ 
ternal to the tile) and global symmetry of the tiling. 

Since the ShellGen algorithm satisfy the decoration rules, it correctly generates symmetric 
shell structures. 


20 




Figure 11: Predicted shell structures for known viruses, (a) Predicted shell structure for 
Tobacco Necrosis Virus using a polyhedra with h=l, k=0. (b) Predicted shell structure for Nudau- 
relia Capensis Virus using a polyhedron with h=2, k=0. (c) Predicted shell structure for Rice 
Dwarf Virus outer shell using h=l, k=3. (d) Predicted structure for the same virus with h=3, k=l. 
All the predicted models have the correct inter-tile (or inter-protein) interfaces/contact in terms of 
geometry, and have the correct global and local topology; except, the one in (c) which has wrong 
chilarity). 

5.2 Reults 

5.2.1 Reproducing Known Shell Structures 

To verify the effectiveness of the tiling and decoration algorithms, we applied it to predict shell 
structures for some viruses for which the structure of the individual building blocks (proteins), as 
well as the entire shell is known. We show some examples here. 

In Figure [T^a), we present the results of modeling the Tobacco Necrosis Virus which has 60 
proteins on its shell. We templated it based on a polyhedra with h=l, k=0 and the resulting 
computationally predicted shell had less that 5A RMSD error with respect to the known shell. 
This error is considered acceptable in molecular biology community. There is no topological errors. 
We had similar success in predicting the structure of Nudaurelia Capensis Virus which requires 240 
proteins on its shell and the template polyhedra was constructed with h=2, k=0 (Figure [T^b)). 
Finally, Figure [TT|(c-d) show the outcome of predicting the shell of Rice Dwarf Virus which has 780 
proteins. The layout for this can be either a polyhedra with h=3, k=l; or h=l, k=3; the latter 
is topologically incorrect if compared to the shell found in nature. Unfortunately, our geometric 
optimization algorithm and scoring model cannot discriminate between the two. 

5.2.2 Assembling Multiple Sized Shells Using the Same Building Block 

In Figure [T^ we show that our algorithm can easily produce shells of different sizes using the same 
building blocks. Here we used the same protein, but decorated tilings of different complexities and 
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Figure 12: Shells of different sizes using the same protein 


reported the highest scoring models for each size. 


6 Conclusion 

We have characterized a new family of polyhedra with regular faces such that it is isotoxal, isohedral, 
and have exactly 2 types of vertices; as well as a dual family which is isogonal, isotoxal and have 
exactly 2 types of regular faces. We have shown that both of polyhedrons of these families generated 
by unfolding a regular polyhedron onto a lattice in a compatible way, thereby allowing the lattice 
vertices, edges and faces to etch out a tiling on the unfolded polyhedron, and finally folding it 
back again. Further, the compatible ways are specified using only a couple of integer parameters. 
We also provided a deterministic and efficient algorithm for generating such polyhedra of any size 
(determined by the two parameters). We have proved that our construction covers all possible 
polyhedron which satisfies the stated properties. When considering the geometric aspects of the 
generated polyhedra, we characterized the cases where the faces may become non-regular, and 
provided solutions for each case. 

Finally, we point out that our class of polyhedron is not similar to the known families like Cata¬ 
lan solids, Johnson solids, or Archimedean solids. Some Catalan solids like the tetrakis hexahedron, 
triakis octahedron, triakis icosahedron, rhombic dodecahedron, rhombic triacontahedron, or pen- 
takis dodecahedron may seem like they satisfy the properties of almost-regular polyhedron, but 
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actually all of them violate either the global of the local symmetry conditions. Also Archimedean 
solids like the truncated cube can be generated by placing the triangles of a tetrahedron on a 
hexagonal lattice such that the corners of each triangle fall on the centers of 3 faces surrounding a 
single face. Many other Archimedian solids are isogonal and isotoxal, but are none of them (not 
even the truncated cube) are duals of of any almost-regular polyhedron. 

The characterization and construction would greatly aid computational analysis and modeling 
of dome or spherical shaped objects with symmetry and regular tiles, which is relevant for anti-viral 
drug design, designing nano-cages for drug delivery and cancer therapy, designing easy to assemble 
masonry structures etc. Additionally, the regular tilings we produce promises to be an interesting 
template for arbitrarily refined meshing of manifolds, by diffeomorphic mappings between the 
manifold and a sphere. 
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